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^ ' or minimization of different physical quantities. We present solutions which 
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Optimal laser control in quantum systems is a problem of fundamental importance for 
atomic, molecular, solid-state and chemical physics which has attracted much attention in 
the last 10 years. Any system exhibiting quantum coherence can be subject to optimal 
control, which basically consists in the manipulation of the quantum dynamics by external 
time-dependent laser fields. A general procedure for this manipulation is described by the 
optimal control theory (OCT) [1,2]. The time-envelope of the external field is optimized by 
experimental or numerical pulse-shaping techniques in order to drive the wave function ip{t) 
of the system to fit a target state 0o at a particular control time t control- Thus, the optimal 
shape Voptit) of the external field achieves the maximization of the quantity | {4'o\4'{'tcx>ntroi))\- 
Since the pioneering work by Hudson and Rabitz [3] , different approaches to optimal control 
have been proposed and experimentally applied in various contexts [1,2,4-10]. The ma- 
jor problem of the current theoretical description of optimal control is that the resulting 
equations are of high complexity and must be solved numerically. Therefore, neither the 
experimental realization nor the theoretical description of optimal control guarantee that 
the result of the optimization corresponds to the true global extremum of the control prob- 
lem considered. As a functional of the pulse shape, the physical quantity to be maximized 
or minimized represents a hyper surface in a multidimensional space, which might exhibit 
many local extrema in which the experimental or numerical procedure can get trapped. 

The only way to extract the global extremum from among the multiplicity of local ex- 
trema is by finding the analytical solution to the control problem. In this paper we present 
for the first time analytical results for the optimal pulse shapes leading to the global extrema 
of different optimal control problems. 

We concentrate on physical situations which can be described by two-level systems. Using 
our approach, we give explicit pulse shapes for inducing the maximization of population 
transfer between two levels and for the achievement of self-induced transparency under the 
constraint of fixed pulse energy. 

We consider a physical quantity Q which is only nonzero when the system is in the 
excited state, and which we wish to maximized or minimized within the control interval 
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[TojTi]. For this purpose an external field of optimal shape has to be applied. Q can be 
seen as a fitness function. We write it as a time average, over the time interval [To,Ti], of 
the form 

Q = r Q{t)dt. (1) 

Jto 

Q{t) is a fitness density, which is a functional of the density matrix p{t) of the system and, 
consequently, also of the external field V{t) cos{u;t). 

Our approach to derive an equation for the optimal field shape is based on a Lagrangian 
of the form [11] 

L= r cdt= r T{t) + iz{t)^p{t)dt + A / ' Q[p{t)]dt +Xi r v\t)dt, (2) 

where C is the Lagrangian density, A and Ai are Lagrange multipliers and r{t) a Lagrange 
multiplier density [12]. 

The first term in Eq. (2) ensures that the density matrix satisfies the quantum Liouville 
equation with the corresponding Liouville operator Z(t). The second term explicitly includes 
the description of the optimal control and refers to a physical quantity to be optimized 
during the control time interval. The third term in Eq. (2) represents a constraint on the 
total energy Eq of the control field [11]. The Lagrangian (2) is the basis of our control theory, 
since it allows the derivation of the equations to be fulfilled by the control field. 

As mentioned before, one of the purposes of this paper is to find analytical solutions for 
the "standard" optimal control problem. This means, we search for the maximization or 
minimization of Q{tcontroi) at a particular control time tcontroi- Note that we can treat this 
problem as a particular case of the theory presented above. We only need to modify the 
fitness density Q by introducing a delta function as follows 

Qipii)]^^ - tcontroddt ^ Q[p{tcontrol)], (3) 

JTq 

where tcontroi lies within the time interval [7o,7i]- 

We now apply the theory described above to a two-level quantum system with energy 
levels ei and 62, interacting with an external control field of the form V{t) cos{u!t), where cu is 



the carrier frequency. If the resonant condition cu — €2 — 61 and adiabaticity criterion \ Vu!\ <^ 
\Vf [13] for the control fields are satisfied, one can use the Rotating Wave Approximation 
(RWA) to derive the Liouville equation for the density matrix 

.dpn 



dt 



jJiV{t){pi2 - P21) + ?7lP22, 



dp22 

i-g^ = pV{t){p2i - P12) - hlP22, (4) 

^ '"^^^^^'^22 - Pll) - «72Pl2, 

where p is the dipole matrix element of the two-level system and 71, 72 are relaxation and 
dephasing constants, respectively. Here we use the notation pi2 — pi2exp{iu!t) and P21 — 
P21 exp{—iu!t) . pii and P22 correspond to the instantaneous occupation of the ground and 
excited state, respectively. Note, that pu + P22 = 1 and p2i = pl2- We set the initial 
conditions as pn = 1, P22 = P12 = P21 = 0. 

As a first application we address the phenomenon known as self-induced transparency 
(SIT) [14,15]. The problem consists in finding a temporal pulse shape for which a light 
pulse entering a material propagates without significant losses. This effect has been studied 
theoretically using different approaches [14,15]. However, it has never been considered so 
far as an optimal control problem. We show below that SIT can be viewed as the search for 
the optimal pulse shape for which losses during propagation are minimized. 

In order to solve the problem from the perspective of optimal control, one can assume, as 
usual, that materials where SIT occurs are collections of inhomogeneously broadened two- 
level systems [14,15]. If we assume, in addition, that the material is optically thin (thus, we 
neglect changes of the field along the spatial axis), we can describe the phenomenon using 
the theory presented above and obtain analytical results. The losses of such a system are 
proportional to the average occupation of the upper level 

^loss OC Q22 = / P22{t)dt. (5) 

Thus, we search for the minimization of the physical quantity Q221 i-e., the integral of the 
occupation of the excited level over the time interval [ro,ri]. Since systems exhibiting SIT 



are characterized by long life-times of the excited levels, we assume 71,27" <^ 1, where T 
is a characteristic time during which the system is in the exited state. T is proportional 
to the inverse Rabi frequency induced by the control field. This assumption permits us to 
neglect relaxation and dephasing effects within the control interval and write the solution 
for occupation of the upper level p22it) as 

P22{t) = sm\e{t)) , (6) 

where the pulse area 9 is defined as 

e{t)^^Ji [ dt'V{t'). (7) 

Thus, the Lagrangian density is a function of the pulse area and its first derivative C{9,9) 
and is given by 

£ = sin2(^(f)) + A^. (8) 
A* 

In Eq. (8), and also in the rest of the paper, we omit the Liouville term since we include 
the analytical expression (6) for the density matrix which solves the Liouville equation. The 
Euler- Lagrange equation derived from the above Lagrangian density reads 

2Xe(t) - iJ^ sm{2e{t)) ^ 0. (9) 

Eq. (9) is of the second order and requires two boundary conditions. We consider for 
simplicity an infinitely large control interval t G (— oo, oo) with natural boundary conditions 
^(— oo) = V{—oo) — V{+oo) — and ^(+00) = tt. This condition imphes that the system 
is excited from and de-excited to the ground-state, remaining there after the interaction 
with the control field. It is clear that these are the boundary conditions compatible with a 
minimization of Q22 (see Fig. 1). Eq. (9) can be integrated analytically (it is mathematically 
equivalent to the pendulum equation). The resulting optimal field envelope Vopt{t) is given 
by the expression 

Voptit) = {^/Xcos^t/^/X/J^)~\ (10) 
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The Lagrange multiplier A is determined from the normalization condition V'^{t)dt — Eq 
for the pulse energy and is given by A = 4/(/i£'o)^. 

Eq. [10] represents the well known soliton solution for the pulse shape and is shown 
in the inset of Fig. 1. This result shows that the sohton wave does not only propagate 
without shape changes [16] but it also minimizes the energy losses, which arc proportional 
to Q22 in the limit of the weak relaxation and dephasing. This fact can be clearly shown in 
Fig. 1, where we show the integrated population p22{t')dt' over the control interval for 
the optimal pulse shape and for a square pulse having the same area and the same energy 
as the optimal. From the figure it is clear that (522, and therefore the losses, is smaller for 
the soliton pulse shape. 

Using other asymptotic values of the pulse area ^(+00) = tt, with N = 2,3, one 
can immediately reproduce 27r, Stt, . . . soliton solutions. These soliton shapes are the optimal 
pulses corresponding to different values of the pulse energy, which are, of course, higher than 
the energy of the shape of Eq. [10]. 

Now we turn to the "standard" optimal control problem at fixed time and show how 
the analytical solution to the problem of maximization of an objective at a particular time 
tcontroi arises naturally as a limiting case of our theory. If the quantity to be maximized is, for 
instance, the population of the upper level at time tcontroi, P22{t control), then the Lagrangian 
density becomes, with the help of Eq. (3) 

JC-S = P22{t)S{t - tcontroi) + XO\t) / , (11) 

where the delta function 6{t — tcontroi) accounts for the modification of the fitness density. 
The optimal pulse shape can only be obtained analytically for 71 = 72 = 0, i.e., if relaxation 
and dephasing effects are neglected. In this case the corresponding Euler- Lagrange equation 

reads 

2X9{t) - nH{t - tcontroi) sm{29{t)) = 0. (12) 

Integrating Eq. (12) one obtains the pulse area as a linear function 6(t) = At + B. By 
substituting the boundary conditions ^(0) = 0, 0{tcontroi) — 7i'/2, we find that the solution of 
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Eq. (12) is a field with a constant amplitude 

Voptit) = , (13) 

^H' ''control 

with energy Eq — tt'^ /{4:iiHcontroi), measured in the time-interval [0, tcontrol]- This result 
reflects the fact that, among all pulses with area equal to 7r/2, that which minimizes the 
energy has time-independent shape. It must be pointed out that this analytical solution 
corresponds to the global extremum of the Lagrangian as long as the RWA is apphcable. 

It is important recall that Eq. (13) is a new result and should not be confused with 
the trivial fact that a 7r/2 pulse, when limited to a constant amplitude, produces complete 
inversion. 

In Fig. 2 we compare the analytical result of Eq. [13] with a numerical solution obtained 
by Zhu et al. using OCT [2]. In Ref. [2], the optimal field to induce population inversion 
between two levels of the Morse potential at a particular time was calculated using an 
iterative numerical technique to integrate the OCT equations. The obtained field consists of 
a single frequency (resonant with the difference of the level-energies) and a time-dependent 
amplitude, shown as a dashed line in Fig. 2. In order to reproduce the same physical 
situation and use the same parameters as in Ref. [2] we calculated the dipole matrix element 
H as n = {ipo\fiipi), where jl{r) = fiore~^^^°,iJ,o = 3.088, and ro = 0.6. ipo{r) and ipi{r) 
correspond to the ground and the first excited state eigenfunctions of the Morse potential, 
which is given by V{r) = Do(exp(-/3(r - r*)) - 1)^ - Dq, with Dq = 0.1994, /3 = 1.189 and 
r* = 1.821 [2]. Thus, we used Eq. [13] to determine the magnitude of the constant optimal 
amplitude. Our analytically calculated optimal field, which corresponds in fact to the true 
global extremum, is shown as a thin sohd line in Fig. 2. Note that the shape obtained by 
Zhu et al. is close to that of the optimal field obtained by us. However, it is clear that 
it does not correspond to the global extremum of the problem. Moreover, the numerically 
determined shape shows a slight asymmetry, which gives rise to a broadening of its Fourier 
spectrum, i.e., to a less effective coupfing to the two- level system. 

This example shows that our analytical approach can be used to check the ability of 
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different numerical methods to avoid local extrema. 

Summarizing, we have obtained analytical solutions for the optimal shape of external 
fields to control populations in two-level systems over finite time-intervals and at particular 
control times. Our obtained optimal shapes constitute the global extrema of the control 
problem, in contrast to previous numerical solutions. Our results can be used as a basis 
to solve optimal control problems in materials which are well described by collections of 
two-level systems. 

This work has been supported by the Deutsche Forschungsgemeinschaft through SFB 
450 and by the Spanish MCyT through BFM2002-03241 and the Program Ramon y Cajal. 
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Time (a.u.) 

FIG. 1. Integrated occupation of the upper level P22{t')dt' for the optimal (soliton) 
pulse (solid line) and for a square pulse (dashed line), shown in the inset figure. Note that 
Q22 = P22(i')'^*' is proportional to the losses Tioss- Inset figure: pulse shape for the opti- 
mal control field Vopt{t) to achieve self-induced transparency through minimization of the losses 
in a two-level system (solid line), and for a square pulse having the same pulse area and energy 
(dashed line). 
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FIG. 2. Thin solid line: analytical solution for the optimal control field Vopt{t) cos {cot) , with 
energy Eq = ir^ / {^n^tcontrol) (see text) to produce inversion of the population at tcontrol = 30000 
(a.u.). The optimal amplitude Vopt{t) is time-independent [see Eq. (13)]. Dashed line: numerical 
result Voct{t) for the field amplitude for the same problem and the same system-parameters obtained 
in Ref. 2 using optimal control theory. 
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